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Abstract 

^ I We present a search for/3 + /EC double beta decay of 120 Te performed with the CUORICINO experiment, an array of Te02 cryogenic 
bolometers. After collecting 0.0573 kg-y of 120 Te, we see no evidence of a signal and therefore set the following limits on the half- 
00 life: T° v /2 > 1 .9 • 10 21 y at 90% C.L. for the Ov mode and T 2 ^ > 7.6 ■ 10 19 y at 90% C.L. for the 2v mode. These results improve the 
'xt" existing limits by almost three orders of magnitude (four in the case of Ov mode). 



C> The discovery of neutrino oscillations JH proved that neu- 
T j~^ trinos are massive, but there are fundamental issues that os- 
cillation experiments cannot address: measuring the absolute 
k> neu tri no mass and determining whether the neutrino is the an- 
• jjj jiparticle of itself, thus being of Majorana nature, or not. To 
C3 answer these questions it is necessary to look for neutrinoless 
double beta, Qv/3/3, decays, which would bring conclusive evi- 
dence of the Majorana nature of the neutrino and whose decay 
rates constrain the absolute neutrino mass Q2Q . 

Double beta decays can occur by either emitting two elec- 
trons or two positrons. In the latter case, either of the 
positron emissions can be replaced by an electron capture (EC). 
While fS~p~ decays have the largest expected rates, /3 + /EC and 
j6 + j6 + decays provide clear signatures from the 51 1-keV annihi- 
lation gamma rays. Energy and momentum conservation in the 
EC/EC decay requires an extra radiative process, reducing the 
rate by several orders of magnitude. 

This paper reports on a search for /3 + /EC decays 120 Te — » 120 
Sn + e + and 120 Te ^ I20 Sn + e + + 2v with the CUORICINO 
experiment, an array of Te02 cryogenic bolometers at the Gran 



Sasso National Laboratories. The 120 Te isotope has been only 
minimally investigated from a theoretical point of view. No cal- 
culations of the half-life of Ov/? + /EC decay of 120 Te are available 
for comparison with our result. The predictions for the same 
decay mechanism in other nuclei, assuming the effective Majo- 
rana mass (m v ) = 1 eV, range between 10 26 -10 27 y JUj. The 
larger value is mostly due to the need to have simultaneously a 
decay and a capture and to the reduced phase space available. 
The only reference for the two neutrino mode 15[] yields a theo- 
retical value for the half-life of 4.4 • 10 26 y. 

In recent years, experimental limits on the 120 Te decays have 
been set using an array of CdZnTe detectors located at the Gran 
Sasso Underground Laboratory |0, 01 and a HPGe detector at 
the Modane Underground Laboratory QS^ l2D- The present best 



limits in the literature are T°^ 2 



> 4.1 • 10" y J7] and T 



1.9-10 17 y[8D. 



,(0v+2v) 
1/2 



The search strategy is the following: the emitted positron 
carries a kinetic energy up to K„ MA = Q - 2m c c 2 - E/,, where E/, 
is the binding energy of the captured electron within the atomic 
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Figure 1 : A sketch of the CUORICINO assembly showing the tower hanging 
from the mixing chamber, the various heat shields and the external shielding. 



shell and Q= (1714.8 ± 1.3) keV H3l is the difference in 120 Te 
and 120 Sn atomic masses. The electron capture is most likely 
to occur from the K shell, whose binding energy is 30.5 keV. 
The ratio of L-capture to K-capture for most elements is around 
10% (12% for 120 Sb -> 120 Sn EC decay) HI. In the following 
we will always assume a capture from the K shell. 

The bolometer where the decay occurs will see the deposition 
of both the binding energy and the kinetic energy of the positron 
(maximum energy: Eq = K max + Eb — Q - 2m e c 2 = 692.8 keV, 
independent of E/,)^ Once at rest, the positron annihilates with 
an electron and two photons with an energy E y = 5 1 1 .0 ke V are 
emitted. These photons can interact with the same bolometer or 
by a nearby one, or they can escape undetected. 

The analysis presented here searches for the signatures with 
the best signal-to-noise ratio. This is not the case for the events 
where the entire energy is deposited inside the detector, due to 
the low detection efficiency (see Table [TJ. For the Ov mode, 
where the positron is monochromatic, this means the coinci- 
dence between a bolometer with an energy deposition consis- 
tent with E y and one with either Eg or Eq + E y and the co- 
incidence of one bolometer with a signal of Eq and two other 
bolometers with a signal of E y . For the 2v mode, where the 
positron is emitted with a continuum of kinetic energies be- 
tween and Kmax, this means the triple coincidence of one 
bolometer with a signal between Et and Eq and two others with 
a signal of E y . 

1. The CUORICINO detector 



capacity 2.3 • 10~ 9 J/K. To measure temperature variations cor- 
responding to few ke V (AT ~ 0. 1 fiK/ke V) heavily doped high- 
resistance germanium thermistors (NTD, Neutron Transmuta- 
tion Doped) are glued to each crystal. 

The CUORICINO detector consists of 62 Te0 2 crystals ar- 
ranged in 13 planes. Each of the upper 10 planes and the lowest 
one consists of four 5x5x5 cm 3 Te02 crystals, while the 1 1th 
and 12th planes have nine, 3x3x6 cm 3 crystals. All crys- 
tals have natural isotopic abundances except four of the smaller 
crystals, two of which are enriched to 82.3% in 128 Te and two 
to 75% in 130 Te. 

The natural abundance of 12() Te is 0.096% HI, so the 39.4 kg 
of the CUORICINO experiment (enriched crystals are not in- 
cluded) contain N m = 1.43 ■ 10 23 nuclei of 120 Te . 

The experiment is shielded with two layers of lead of 10 cm 
minimum thickness each. The outer layer is made of common 
low radioactivity lead, while the inner layer is made of special 
lead with a low activity of 21() Pb. The electrolytic copper of the 
refrigerator thermal shields provides an additional shield with a 
minimum thickness of 2 cm. An external 10 cm layer of borated 
polyethylene was installed to reduce the background due to en- 
vironmental neutrons. The detector itself is shielded against the 
intrinsic radioactive contamination of the dilution unit materials 
by an internal layer of 10 cm of Roman lead 11411 . located in- 
side the cryostat immediately above the tower. The background 
from the activity in the lateral thermal shields of the dilution re- 
frigerator is reduced by a lateral internal 1 .4 cm thick shield of 
Roman lead. Another 8 cm lead shield is located at the bottom 
of the tower. The refrigerator is surrounded by a Plexiglas anti- 
radon box flushed with clean N2 from a liquid nitrogen evapo- 
rator and is also enclosed in a Faraday cage to eliminate elec- 
tromagnetic interference. A sketch of the assembly is shown in 
Fig. ID 



Signature [energies in keV] 


ft 


£[%] 


(30.5-692.8) 


1 


3.00 ± 0.02 


(30.5-692.8) + 511 


2 


3.40 ± 0.02 


(30.5 -692.8) + 511 + 511 


3 


0.45 ± 0.01 


(541.5- 1203.8) 


1 


16.28 + 0.04 


(541.5 - 1203.8) + 511 


2 


6.23 ± 0.03 


(1052.5- 1714.8) 


1 


10.04 ± 0.03 



Table 1: Signatures of I20 Te /3 + /EC decay in an array of TeC>2 detectors and 
their corresponding multiplicity (/i), that is the number of detectors with an 
energy deposition above threshold. The detection efficiency for the Ov mode 
in CUORICINO is reported in the last column (e). We denote with the + sign 
the coincidence of energies released in different detectors. For the Ov mode the 
energy released in the detector where the decay occurred corresponds to the 
upper bound of the interval. The errors are statistical only. 



The CUORICINO experiment is detailed in Ref. 112]. 
Briefly, it is an array of Te02 crystals acting as cryogenic 
bolometers at a working temperature of 8-10 mK and with heat 



Here and in the following we assume that the X-rays following the EC do 
not escape from the crystal where the decay occurs. 



For the present analysis, the full CUORICINO statistics (data 
collected between May 2004 and May 2008) for a total expo- 
sure of 0.0573 kg-y of 120 Te is used. The total energy spectrum 
of all detectors is shown in Fig. [2] Several peaks of radioactive 
isotopes are clearly visible, the most prominent are labeled. 
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Figure 2: Total energy spectrum of all CUORICINO detectors. The most prominent peaks are labeled and come from known radioactive sources such as: e + e 
annihilation (1), 214 Bi (2), 40 K (3), 208 T1 (4), 60 Co (5) and 228 Ac (6). 
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Figure 3: Resolution fits at E=51 1 keV (left) and E=1238 keV(right). The fit function is a Gaussian plus linear background. 



2. Data acquisition and analysis 

CUORICINO data are divided in datasets, each one being a 
collection of about a month of daily measurements. Routine 
calibrations are performed at the beginning and at the end of 
each dataset using two wires of thoriated tungsten inserted in- 
side the external lead shield. The signals coming from each 
bolometer are amplified and filtered with a six-pole Bessel low- 
pass filter and fed to a 16-bit ADC. The signal is digitized with 
a sampling time of 8 ms. With each triggered pulse, a set of 5 12 
samples is recorded to disk. The typical bandwidth is approxi- 
mately 10 Hz, with signal rise and decay times of order 30 and 
500 ms, respectively. More details of the design and features of 
the electronics system are found in Ref. 1 15]. 

Each bolometer has a different trigger threshold, optimized 
according to the bolometer's typical noise and pulse shape. The 
trigger rate is time and channel dependent, with a mean value 
of about 1 mHz. The amplitude of the pulses is estimated by 
means of an Optimal Filter technique 11611 . The gain of each 
bolometer is monitored by means of a Si resistor of 50-100 k£2 
attached to it that acts as a heater. Heat pulses are periodically 
supplied by an ultra-stable pulser [17] that sends a calibrated 



voltage pulse to the Si resistors. Their Joule dissipation pro- 
duces heat pulses in the crystal with a shape which is almost 
identical to the calibration y-rays. 

3. Search strategy 

As detailed in the introduction, a /3 + /EC decay of 120 Te in 
an array of Te02 detectors releases an energy up to Eq = 
692.8 keV in the bolometer where the decay occurs and, if the 
annihilation gammas do not escape undetected, one or two ad- 
ditional energy deposits of E y = 5 1 1 .0 ke V in either the same 
or a nearby bolometer. There are therefore several distinctive 
signatures as listed in Tab. Q] 

For the Ov mode the energy released in the detector where 
the decay occurred corresponds to the upper bound of the inter- 
val. The detection efficiencies for the Ov mode were estimated 
by means of a GEANT4 simulation of the CUORICINO setup 
(see last column of Tab. [TJ where the decays are located uni- 
formly within all non-enriched detectors and the decay products 
are emitted isotropically. We always assume that the binding 
energy of the captured electron is released within the detector 
where the decay occurs. The efficiency estimate includes the 
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Figure 4: Energy spectra of coincidences with a single 511 keV event (left) and of coincidences with two 511 keV photons (right). The most prominent peaks are 
labeled and come from the single escape and the double escape of known lines from the following isotopes: 208 T1 (1), 214 Bi (2) and 40 K (3). In the spectrum of 
double coincidences (left) the 511 keV peak is also visible (4). 



dead time evaluated separately for each detector. To account 
for the different energy resolution of CUORICINO detectors 
we assigned to each detector in the simulation a FWHM given 
by the weighted mean of the FWHMs calculated in all CUORI- 
CINO calibration runs. The highest detection efficiencies cor- 
respond to the cases where one or both of the electron-positron 
annihilation photons are fully absorbed in the detector where 
the p + /EC decay occurs. This is consistent with the relatively 
large size of the CUORICINO detectors and the mean free path 
of a 5 1 1 ke V photon in Te02 (1.9 cm). These signatures are cal- 
culated to account for nearly half of all decays. The remainder 
involve only partial energy deposition of the 511 keV gamma 
ray energy and therefore are more difficult to distinguish from 
background. Due to the extremely short range of positrons 
in CUORICINO crystals, the efficiencies estimated for the Ov 
mode are also valid for the 2v mode, apart from small correc- 
tions (see below). Considering also the expected background, 
we limited our analysis to the signatures that feature a coinci- 
dence of two or three events for the Ov mode and of three events 
for the 2v mode. 

We consequently search only for double or triple coinci- 
dences where one or two of the energy depositions are in the in- 
terval +2.5<x of E y = 511 keV, where cr — 1.7 keV, as estimated 
in the inclusive energy spectrum (see Fig. [3j. The coincidence 
window is 100 ms. The probability of accidental coincidence 
is estimated from the measurement of the single crystal rate (~ 
1.6 mHz) to be around 0.7% for double coincidences between 
detectors in the same plane or in adjacent planes. For triple 
coincidences it is even less. 

The energy spectra of the events in coincidence with one or 
two 511 keV photons are shown in Fig. [4] The structures in 
both spectra have been identified and correspond to single or 
double escape lines from known radioactive lines present in 
the CUORICINO background spectrum. The continuum back- 
ground is due both to accidental coincidences and true coinci- 
dences in which the energy deposition of the event that accom- 
panies the 511 keV gamma(s) is not complete. 



4. Results of the 0vj8 + /EC decay search 

In the spectra of Fig. |4] we search for a signal with mean 
energy 692.8 keV (spectra of double and triple coincidences) 
or 1203.8 keV (spectrum of double coincidences only). The 
expected resolution at 1203.8 keV is estimated on the 214 Bi line 
at 1238 keV and found to be cr = 1 .67 + 0.05 keV (see Fig.|3), in 
agreement with the value of 1 .7 keV found for the 511 keV line 
and consistent with being constant over the energy region of 
interest. 

Fig. [3] shows the energy regions of interest of the double 
and triple coincidences spectra. The observed lines, as ex- 
plained above, correspond to single- and double-escape peaks. 
Specifically, the measured spectra show the single escape lines 
of the 1173.2 keV transition from 60 Co at 662.2 keV and of 
the 1729.6 keV transition from 2I4 Bi at 1218.6 keV plus the 
double escape lines of the 2204.06 keV transition from 2I4 Bi 
at 1182.06 keV, of the 1764.5 keV transition from 214 Bi at 
742.5 keV and again of the 1729.6 keV transition from 2I4 Bi 
at 707.6 keV. The observed SE and DE peaks produce a neg- 
ligible contribution to the background in the energy regions of 
interest. 

The search for the 0vy6 + /EC decay of 120 Te on the spectra in 
Fig. [5] is performed with unbinned maximum likelihood fits. In 
the likelihood fit the background is therefore parametrized with 
the sum of a flat component and a Gaussian for each of the ex- 
pected escape lines, with mean values fixed to the correspond- 
ing known energies and the width fixed to cr — 1 .7 keV. The sig- 
nal is also parametrized with a Gaussian with mean value fixed 
to the expectations, that is 692.8 keV or 1203.8 keV, depending 
on whether or not one of the annihilation photons is absorbed in 
the same crystal where it is emitted. The width of the signal line 
is fixed to 2.2 keV which is obtained by summing in quadrature 
the sigma (1.7 keV) , the error on the Q-value (1.3 keV), and 
the uncertainty on the energy scale (0.4 keV) as estimated from 
the calibration fij^. In summary for each spectrum the fitted 



2 From a bayesian point of view (i.e. our approach), you should sum the 
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Figure 5: Energy spectra, in the regions of interest, of coincidences with a single 
511 keV event [two instances] (top) and of coincidences with two 511 keV 
photons (bottom). The arrows indicate the energy of the expected signal. Fit 
results, as explained in the text, are overlaid. 



quantities are the number of events in the signal, the flat back- 
ground and the escape lines. The curves resulting from the fits 
are overlaid in Fig.[5]and the corresponding fitted quantities are 
reported in Tab. [2] 

The upper limit on the number of observed signal events is 
extracted by means of a Bayesian procedure. By integrating 
with a Monte Carlo technique the likelihood over the nuisance 
parameters, we calculate, for each signature, the posterior prob- 
ability density function (pdf) for the parameter N s i g and define 
our 90% C.L. upper limit as the value where its integral reaches 
the 90% of the total area (see Tab.©. 

To combine the results coming from the three signatures, 



likelihoods for different Q- values weighting them by the probability of that Q- 
value being correct. This is equivalent to convoluting the peak position with a 
gaussian, and therefore to a gaussian whose width is the sum in quadrature of 
the uncertainty on the Q- value and the resolution. 



we calculate with the same procedure the pdfs for the param- 
eter N S ig/e tot , the number of signal events divided by the cor- 
responding efficiency. The efficiency is calculated as e tot = 
£ x (CnoisetheatT \ where € is the efficiency tabulated in Tab. Q] 
[i — 2, 3 is the multiplicity of the given signature from the same 
table, e no i se =99.1% accounts for the loss of signal due to noise 
as estimated on heater pulses, and Cheat - 97.7% accounts for 
the dead time induced by the presence of heater. This proce- 
dure, together with the inclusion of the uncertainty on the en- 
ergy scale in the fits of the spectra, accounts for the systematic 
errors, which have a negligible impact on the result. 

From the combined pdf, we extract a 90% C.L. upper limit 
on the total number of double beta decays regardless of their 
detection, n B = 100. We can then set a limit on the half-life of 
0v/3 + /EC decay of 12l) Te 



nOv 

1/2 



>ln2Nm— = 1.9- 10 21 y 



n B 



(1) 



where Npp is the number of Te nuclei and T is the live time. 

5. Results of the 2v/? + /EC decay search 

The energy spectrum of the events in coincidence with two 
511 keV photons in the region below 692.8 keV contains 15 
events (see Fig. [§}. We can exclude from our analysis the re- 
gions +3cr around the energies where we expect the DE peaks 
from known radioactive gamma lines in the CUORICINO in- 
clusive spectrum. These energies are indicated with red arrows 
in Fig. [6] and listed in the caption. We have listed only the 
gamma lines that are expected to contribute with at least one 
event. This expectation is based on the comparison with the 
40 K line that produces 5 events at 438.8 keV in the experimen- 
tal spectrum of triple coincidences (see Fig. [6] line labeled as 
1). After the subtraction, only 8 events remain in the spectrum 
from threshold to end-point. These events can be due to acci- 
dental coincidences or to true coincidences (double escape of 
already compton scattered gammas). To estimate the first com- 
ponent we looked at the spectra of events in triple coincidence 
with the side-band of 5 1 1 keV (left side-band: from 470 keV 
to 502.5 keV, right side-band: from 519.5 keV to 560 keV) cor- 
rectly normalized to the experimental spectrum of Fig. [6] This 
accounts for 4.3 + 0.5 events. The other component of the back- 
ground cannot be reliably estimated due to limited statistics and 
we cannot discriminate against it. Therefore, we set an upper 
limit assuming conservatively that the remaining events may be 
signal. 

The upper limit on the number of observed signal events is 
extracted by means of a Bayesian procedure. From the likeli- 
hood function for Poisson distributed data with unknown mean 
s and known background b and using as prior pdf a flat prior 
different from zero only for positive values of s, we extract our 
90% C.L. upper limit on the number of observed events as the 
value where the integral of the posterior pdf reaches the 90% of 
the total area 1180, n B = 9.04. We can then set a limit on the 
half-life of 2v/37EC decay of I20 Te 

' (2) 



r 2 ; 2 >in2JV 



n B 



0.76- 10 2 % 



spectrum 




Nfiat 


N L i 


N L2 


N H 


90% U.L. 


I 


1.7+3.5 


214+16 


17.5+5.3 


14±5 


34±6 


8.5 


II 


0.3+3.1 


78±10 


26+6 


N/A 


9.5+4.3 


7.25 


III 


0.0±0.5 


1 + 1 


N/A 


N/A 


8±3 


2.62 



Table 2: Results of the unbinned likelihood fits to the energy spectrum of the double coincidences around E = 692.8 keV ("I", shown in the top plot of Fig. [5), 
around E = 1209.8 keV ("II", shown in the middle plot of Fig. [5), and the triple coincidences around E = 692.8 keV ("III", shown in the bottom plot of Fig. [3). 
Here N s j g estimates the number of signal events, Nfua the total number of flat background events in the fit interval, and N12 are me number of events in the 
escape peaks of lower energies and Nh is the number of events in the escape peak of higher energies. Finally, the 90% C.L. upper limits on the number of observed 
events are reported. 




Figure 6: Energy spectrum of triple coincidences with two 511 keV photons in 
the energy region where we expect a signal from the 2vj8 + /EC. The red arrows 
indicate the energies where we expect the DE peaks from known radioactive 
gamma lines in the CUORICINO inclusive spectrum. For each of these energies 
we subtract an energy region of ±3cr. The lines that can give a double-escape 
peak between 50 keV and 692.8 keV are 1120.3 keV and 1238.1 keV of 214 Bi 
(3), 1173.2 keV and 1332.5 keV of 60 Co (2) and 1460.8 keV of 40 K (1). 



where Npp is the number of 120 Te nuclei and T is the live time. 
The efficiency is calculated as e tot = fX(e„«e), eo( ) 3 x(l -e corr ), 
where e, e no ; se and ei, ea , have been defined above and e corr = 
12.5 + 0.5% is the correction to be applied to the detection effi- 
ciency of the Ov mode and accounts for the fact that we are not 
sensitive to the portion of the positron spectrum between 30.5 
and 50 keV (threshold) and that we subtract five 10 keV energy 
regions from the experimental spectrum, as explained above. 
Since we do not know the exact shape of the positron spec- 
trum from the 2v/J + /EC decay of 120 Te the estimation of e co „ 
has been performed on a standard /3 spectrum with end-point at 
692.8 keV and the error is estimated from the comparison of 
e corr estimated as explained above and e corr estimated on a flat 
spectrum from to end point. Systematic errors on the efficien- 
cies have a negligible impact on the result. 



with respect to the constraints from previous experiments. The 
limits obtained with CUORICINO could be improved in the fu- 
ture by an additional ~2 orders of magnitude with the CUORE 
experiment lfl9ll because of increased mass, higher coincident 
detection efficiency, and lower backgrounds. 
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6. Conclusions 

We searched for double beta decays /i + /EC of I20 Te in the 
Te0 2 cryogenic bolometers of CUORICINO, using 0.0573 kg-y 
of 120 Te of data. We see no evidence of a signal and therefore set 
new limits on the half-life for Ov and 2v decay T ( ^ 2 > 1 .9 • 10 21 y 
and T 2 J 2 > 0.76 • 10 20 y extending the exclusion region by al- 
most three orders of magnitude (four in the case of Ov mode) 
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